Condensate Fluctuations in Trapped Bose Gases: 
Canonical vs. Microcanonical Ensemble * 



Martin Holthaus and Eva Kalinowski 

Fachbereich Physik der Philipps-Universitdt, Renthof 6, D-35032 Marburg, Germany 



Klaus Kirsten 

Universitdt Leipzig, Institut fiir Theoretische Physik, Augustusplatz 10, D-04109 Leipzig, 

Germany 
(April 15, 1998) 



Abstract 

We study the fluctuation of the number of particles in ideal Bose-Einstein 
condensates, both within the canonical and the microcanonical ensemble. Em- 
ploying the Mellin-Barnes transformation, we derive simple expressions that 
link the canonical number of condensate particles, its fluctuation, and the 
difference between canonical and microcanonical fluctuations to the poles of 
a Zeta function that is determined by the excited single-particle levels of the 
trapping potential. For the particular examples of one- and three-dimensional 
harmonic traps we explore the microcanonical statistics in detail, with the 
help of the saddle-point method. Emphasizing the close connection between 
the partition theory of integer numbers and the statistical mechanics of ideal 
Bosons in one-dimensional harmonic traps, and utilizing thermodynamical ar- 
guments, we also derive an accurate formula for the fluctuation of the number 
of summands that occur when a large integer is partitioned. 
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There is, essentially, only one problem in 
statistical thermodynamics: the distribu- 
tion of a given amount of energy E over 
N identical systems. 

Erwin Schrodinger |I| 



I. INTRODUCTION 

The ideal Bose gas is customarily treated in the grand canonical ensemble, since the 
evaluation of the canonical partition sum is impeded by the constraint that the total particle 
number N be fixed. In contrast, after introducing a variable that is conjugate to A^, the 
fugacity z, the computation of the ensuing grand canonical partition function 'E[z, (3) requires 
merely the summation of geometric series, and all thermodynamic properties of the Bose gas 
are then obtained by taking suitable derivatives of with respect to z or the inverse 

temperature [3. 

There is, however, one serious failure of the grand canonical ensemble. Grand canonical 
statistics predicts that the mean-square fiuctuation {6'^n,^)gc of the z/-th single-particle level's 
occupation equals {nu)gc {{nu)gc + !)• Applied to the ground state 1^ = 0, this gives 

{S^no)gc = (no) gc {{no) gc + I) 

even when the temperature T approaches zero, so that all particles condense into the 
ground state. But the implication of huge fluctuations, {S^no)gc = N{N + 1), is clearly 
unacceptable; when all particles occupy the ground state, the fluctuation has to die out. 

This grand canonical fluctuation catastrophe has been discussed by generations of physi- 
cists, and possible remedies have been suggested within the canonical framework [0-^. Af- 
ter the recent realization of Bose-Einstein condensates of weakly interacting gases of alkali 
atoms |]5|-^ had brought the problem back into the focus of interest [pUj-P^ , signiflcant steps 
towards its general solution could be made. Ideal condensate fluctations have been com- 
puted for certain classes of single-particle spectra, i.e., for certain trap types, both within the 
canonical ensemble , where the gas is still exchanging energy with some hypothetical 



heat bath, and within the more appropriate microcanonical ensemble |10.T^,T^, where it 



is completely isolated. Interestingly, canonical and microcanonical fluctuations have been 
found to agree in the large- A^-limit for one-dimensional harmonic trapping potentials |lT0| , p!3 
but to differ in the case of three-dimensional isotropic harmonic traps |TT||T^. 



Yet, in the true spirit of theoretical physics one would clearly like to have more than 
merely some formulas for condensate fluctuations in particular traps. Can't one extract a 
common feature that underlies those formulas, such that simply inspecting that very feature 
allows one to determine, without any actual calculation, the temperature dependence of 
the condensate fluctuation, and to decide whether or not canonical and microcanonical 
fluctuations are asymptotically equal? 

It is such a reflned understanding that we aim at in the present work. As will be shown, 
the feature imagined above actually exists: It is the rightmost pole, in the complex t-plane, 
of the Zeta function 
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that is furnished by the system's single-particle energies Sy. In order to substantiate this 
statement, we will proceed as follows: We start in the next section by deriving simple 
expressions that relate the canonical number of ground state particles (no)cn; and its mean- 
square fluctuation ((5^no)cn, to Z{P,t). The key point exploited there is the approximate 
equivalence of a trapped Bose gas in the condensate regime to a system of Boltzmannian 
harmonic oscillators. This equivalence, which holds irrespective of the particular form of 
the trapping potential, implies that both {nQ)cn and {5'^no)cn can be expressed in terms of 
harmonic oscillator sums, which explain the emergence of the spectral Zeta function Z{(3,t) 
and can be computed with the help of well-established techniques [0. In Section |T| we 



then evaluate these general canonical formulas for rf-dimensional isotropic harmonic traps, 
where Z{f],t) reduces to ordinary Riemann Zeta functions, and for anisotropic harmonic 
traps, where it leads to Zeta functions of the Barnes type. 

In Section |IV| we compare the canonical statistics of harmonically trapped gases for d = 1 
and d = 3 to their microcanonical counterparts. The strategy adopted there — the calcula- 
tion of microcanonical moments from the easily accessible corresponding canonical moments 
by means of saddle-point inversions — is technically rather cumbersome and certainly not to 
be recommended if one merely wishes to obtain the microcanonical condensate fluctuations 
('5^?T-o)mc; but it explains in precise detail just how the difference between the canonical and 
the microcanonical ensemble comes into play, and why (5^?T,o)mc becomes asymptotically 
equal to (5^no)cn in some cases, but not in others. A convenient expression for the immedi- 
ate determination of (5^no)cn — {S'^no)mc, based again on the spectral Zeta function Z{P,t), 
is then derived in Section 0. The final Section ^ contains a concluding discussion; three 
appendices offer technical details. 

Since we restrict ourselves to non-interacting Bose gases, the main value of the present 
work lies on the conceptual side — after all, the ideal Bose gas ought to be properly un- 
derstood — , but it may well turn out to be of more than purely academical importance: 
After it has been demonstrated now that the s-wave scattering length in optically confined 
condensates can be tuned by varying an external magnetic field |TB|, the creation of almost 



ideal Bose-Einstein condensates might become feasible. If it did, then also the experimental 
investigation of the basic statistical questions studied here, such as the connection between 
the temperature dependence of the condensate fiuctuation and the properties of the trap 
potential, should not remain out of reach. 

Finally, there is still another appealing side- aspect: Since the microcanonical statistics 
of ideal Bosons in one- dimensional harmonic traps can be mapped to the partition theory 
of integer numbers, a natural by-product of our work is a fairly accurate formula for the 
fiuctuation of the number of integer parts into which a large integer may by decomposed. The 
derivation of that formula in Section |V| is a beautiful example for the deep-rooted connection 
between partition theory and statistical mechanics, quite in the sense of Schrodinger's remark 
quoted above. 
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II. CANONICAL DESCRIPTION OF IDEAL BOSE-EINSTEIN CONDENSATES 



We consider a gas of non- interacting Bose particles confined in a trap witli discrete single- 
particle energies Ei, {u = 0, 1, 2, . . .) and stipulate that the ground state energy be equal to 
zero, £o = 0. Starting from the grand canonical partition sum 

oo -I 



i/=0 



where /? = l/{kBT) is the inverse temperature, we have the familiar expansion 

oo 

E{z,P) = ^;.^^e"'^^l](E|iV), (2) 

N=0 E 

with coefficients fl{E\N) denoting the number of microstates accessible to an iV-particle 
gas with total excitation energy E. Combinatorically speaking, fl{E\N) is the number of 
possibilities for sharing the energy E among up to N particles — the number A'ex of particles 
that are actually excited and thus carry a part of E remains unspecified. 

The clear distinction between and A^^cx is the starting point for studying statistical 
properties of Bose-Einstein condensates in gases with fixed particle number. When A'ex out 
of N Bose particles are excited, there remain N — Ngx particles forming the condensate, and 
the corresponding number of microstates (that is, the number of possibilities for distributing 
the excitation energy E over exactly A^ex particles) is given by 

n{E\N,,) - n(E|iVex - 1) = $(iVex|^) . (3) 

Within the canonical ensemble, i.e., if the A^-particle gas is in contact with some heat bath 
of temperature T, the probability for finding Nex excited particles can then be written as 

The expectation value (A'^ex)cn with respect to this distribution yields the canonical ground 
state occupation number, 

{no)cn = iV - {N^x)cn ; (5) 

the canonical mean-square fluctuation of the number of condensate particles is identical to 
the fluctuation of the number of excited particles, 

{S^nQ)cn = (^^iVcx)cn 

= (Ocn - (iVex)L • (6) 

In order to calculate these cumulants, we consider the function 

{l-z)E{z,P) ^ E,^{z,f3), (7) 

which satisfies the identities 
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oo 

' — 'f^m [ ^ ^ , ^ 

N=0 ' ' E 

oo 



oo 



= E E e""'' [^(^1^) - ^(^1^ - 1)] (8) 

N=0 E 

with VL{E\ — 1) = 0. Hence, replacing the summation index by Ai'ex, one finds 

oo 

Sex(^,/5) = E ^''^''Ee-^'''^(^cx|i?). (9) 
Afcx=0 E 

On the other hand we have 

oo 1 

""'-■^'-n ,_,,,,(_,,^) . (10) 

where, in contrast to Eq. (P, the product runs only over the excited states v > 1: The 
grand canonical partition sum of a fictituous Bose gas which emerges from the actual gas 
by removing the single-particle ground state is the generating function for $(A^ex|-E')- Differ- 
entiating this generating function k times, and then setting z = 1, one gets the canonical 
moments [] 



Y.e-'^'i E N^.HNe.\E)] ^ M,{P). (11) 

,Afox=0 



In the customary grand canonical framework, z is identified with the fugacity and linked to 
the ground state occupation number {no)gc by 2; = (1 + l/{no)gc)~^- In that case z remains 
strictly less than unity, thus preventing the ground state factor in Eq. from diverging. In 
contrast, z is no more than a formal parameter in the present analysis, entirely unrelated to 
the ground state occupation; and since the ground state factor is absent in the generating 
function Sex(2,/9), there is nothing to prevent us from fixing z = 1. 

If the sum over iVgx in Eq. ([TT|) did not range from to cxd, but instead from to 
the actual particle number A^, as it does in the canonical distribution (^, then the ratio 
Mi{P)/Mq{P) would be exactly equal to the canonical expectation value {Ncx)cn- But even 
if we do not have an exact equality here, the difference between these two quantities must be 
negligible if there is a condensate. This statement requires no proof, it is a mere tautology: 
a condensate can only be present if those microstates where the energy E is spread over 
all particles are statistically negligible, so that also the microstates that would become 
available if additional zero-energy (ground state) particles were added to the gas cannot 
make themselves felt. Hence, in the presence of a Bose-Einstein condensate we have, for 
small fc, 



^This corresponds to a trick suggested by D. H. Lehmer and taken up by L. B. Richmond for 
calculating A:-th moments of partitions of integer numbers; see Ref. |19|. Note that the asymptotic 
moment formula derived in that paper is not correct for k > 2; the corrected formula is stated in 
Appendix |C|. 
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N 



Y: N^^<I>{N,^\E) = Y: N^.<I>{N,^\E) (12) 

iVex=0 iVex=0 



at least to a very good approximation, which gives both 

Mm 

MoiP) 



- m (13) 



and 
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The approximation ([T2|) , expressing the replacement of the actual condensate of iV — 
{Nex)cn particles by a condensate consisting of infinitely many particles, can be interpreted 
in two different ways. On the one hand, the infinitely many ground state particles may be 
regarded as forming a particle reservoir for the excited-states subsystem. Such an approach 
to computing canonical condensate fluctuations had been suggested as early as 1956 by 
Fierz B]; it corresponds to the "Maxwell's Demon Ensemble" recently put forward by Navez 



et al. 15 



The second interpretation rests on the observation that for /c = the approximation (p!2D 
takes the form 

oo 

Y $(iVex|^) = niE\N) , (15) 

Arex=0 

SO that Eq. (|ll|) gives 

oo -1 

^e-'^^nim)^n^—^^.ZW. (16) 

Since each factor 1/[1 — exp(— /J^j,)] corresponds to a geometric series, i.e., to the canonical 
partition function of a simple harmonic oscillator with frequency e^/h, Eq. ([T6|) states that 
if there is a condensate, so that Eq. (|T^) holds, then the canonical partition function of an 
ideal Bose gas with arbitrary single-particle energies is well approximated by the canonical 
partition function of a system of distinguishable harmonic oscillators, each excited single- 
particle level El, corresponding to an oscillator with frequency Su/h. Thus, for temperatures 
below the onset of Bose-Einstein condensation the thermodynamics of the actual Bose gas 
practically coincides with the thermodynamics of a Boltzmannian harmonic oscillator sys- 
tem, regardless of the specific form of the trapping potential. For this reason, we will refer 
to the approximation ( [T^ ) as the oscillator approximation. For the particular case of a 
three-dimensional isotropic harmonic trapping potential, the quality of this approximation 
has been confirmed in Ref. [^] by comparing the entropy of the actual Bose gas with that 
of its Boltzmannian substitute. 

Within this oscillator approximation, the determination of the number {Nex)cn of excited 
particles, and of the canonical mean-square condensate fluctuation {5'^no)cn = (5^^ex)cn, 
becomes remarkably simple. Doing the derivatives demanded by Eq. (0), we find 
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Mo(/5) = Z{P) (17) 
Mm = Z{P)S,{P) (18) 
M2(/3) = Z(/3)[52(/5) + 52(/3)l , (19) 



with Z{l3) as given by Eq. (|T6D, and 



y=i exp(/3e^) - 1 

oo oo 

EEexp[-/?e.(r + l)], (20) 



i/=l r=0 
oo 



^'^^^ exp(/3e,) - 1 iexp(/5e.) - 1 + \ 



oo oo 



= E E ^ exp[-/?£,r] . (21) 

v=l r=l 

Computing the ratios Mi{(3) / Mq{(3) and M2{(3) / Mq{(3) according to Eqs. (|13D and (|14D, the 
oscillator partition function Z{[3) drops out (and, hence, does not even have to be evaluated 
here!), and we arrive at the appealing relations 

(iVex),„ = Si{f3) (22) 
(5'iV,,),„ = SM . (23) 

Since l/[exp(/9ei,) — 1] = {n^)gc is just the grand canonical expectation value for the occu- 
pation of the z/-th excited state in a partially condensed Bose gas (the fugacity of which is 
z = 1), and {n^)gc {{nu)gc + 1) = (5^^^i/)gc is the corresponding grand canonical mean-square 
fluctuation, the respresentations (|20| ) and (^I]) reveal that the canonical expectation value 
of the number of excited particles equals the grand canonical one, and that the canonical 
mean-square fluctuation of the ground state occupation number can simply be computed 
by adding the grand canonical fluctuations of the excited levels, subject to only the oscil- 
lator approximation. This is precisely what had been anticipated by Fierz [^, and what 
has also been exploited in a heuristic manner by Politzer when investigating the three- 
dimensional harmonic trap. 

For evaluating the sums Si{(3) and S2{P) we employ the Mellin-Barnes integral repre- 
sentation [ p!7| 

= — / dta-'T{t), (24) 
valid for real r > and complex a with Ih (a) > 0. This leads to 

(iVex)cn = E E ^ X__ [^^^(^ ^ 1)], 
1 /-T+ioo oo oo -p/.x 

^ dtEE 



27ri Jr-i^ ^ [(3eu{r + 1)]* 

/ dtr(t)Z(/3,t)C(t), (25) 

J T—ioo 



2Txi Jt- 
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where ({t) = Yl^^iV * denotes the Riemann Zeta function, and we have introduced the 
spectral Zeta function 

oo -| 

that embodies the necessary information about the trap spectrum. In the same way we find 
the remarkably similar-looking equation 

(5'iVex)cn = — / dtT{t)Z{(3,t)C{t-l) . (27) 

It should be noted that interchanging summations and integration in Eq. (^), and in the 
analogous derivation of the canonical fluctuation formula (pTj), requires the absolute con- 
vergence of the emerging sums. Therefore, the real number r has to be chosen such that 
the path of integration up the complex t-plane lies to the right of the poles of both Zeta 
functions. 

Now the temperature dependence of {Nex)cn or {S^Nex)cn is determined by the pole of 



the integrand (p5|) or (|27| ) that lies farthest to the right. Since T(t) has poles merely at 
t = 0, —1, —2, . . . , the decisive pole is provided either by the Riemann Zeta function ({t) or 
({t — 1), respectively, or by its spectral opponent Z{P,t), which depends on the particular 
trap under study |]16|. This competition will be discussed in detail in the following section, 
focussing on harmonic trapping potentials. 



III. ISOTROPIC AND ANISOTROPIC HARMONIC TRAPS 



The evaluation of the canonical relations (25) and (27) reduces to a mere formality if 
the pole structure of the spectral Zeta function ( [2B| ) is known. The simplest examples are 
provided by c?- dimensional isotropic harmonic traps, since then Z{f3,t) becomes a sum of 
Riemannian Zeta functions. Namely, denoting the angular frequency of such a trap by u, 
the degree of degeneracy g^, of a single-particle state with excitation energy uhcu is 

V ' ^^^^ 

so that Z{(3, t) acquires the form 

ZiP,t) = (fShur'Y.^^ (29) 




giving in explicit terms 



Z{P,t) = {I3hu)-X{t) for d = 1 , 

Z(/3, t) = Iphiu)-' [C(t - 1) + C(^)] for = 2 , (30) 

Z{P, t) = iPhcu)-' [C(t - 2)/2 + 3C(t - l)/2 + cm for = 3 . 

We now aim at the temperature dependence of (A^ex)cn and (5^A^ex)cn for temperatures below 
the onset of a "macroscopic" ground state occupation, so that the oscillator approximation 
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retains its validity, but well above the level spacing temperature hu/kB, so that phu <ti 1. 
Such a temperature interval exists if the particle number N is sufficiently large, since the con- 
densation temperature generally increases with N. The desired asymptotic T-dependence 
can then directly be read off from the residue of the rightmost pole of the respective inte- 
grand (|25|) or p7|). Since ({z) has merely a single pole at z = 1, simple and with residue 
+1, namely ^1 



z 



1 



+ 7 



(31) 



for z close to 1, the calculation of that residue is particularly easy if the rightmost pole 
in ( P5| ) or ( P7D is simple. In the case of a double pole we also need the identity 



r'(n) = r(n) il){n) 

= f-7 + E -) 



(32) 



for the Psi function at integer arguments, with 7 ^ 0.5772 denoting Euler's constant. This 
is the only technical knowledge required for computing the number of excited particles, and 
its fluctuation, in a c?- dimensional isotropic harmonic trap within the canonical ensemble: 

(i) For d = 1, the number of excited particles is governed by the double pole at t = 1 which 
emerges since t) now is proportional to whereas the mean-square fluctuation 
is dominated by the simple pole oi C,{t — 1) aX t = 2: 




(33) 
(34) 



(ii) For d = 2, the rightmost pole of Z{(3, t) has moved to t = 2 and thus determines {Nex)cn 
all by itself, but now the product Z{j3,t)((t — 1) provides a double pole that governs 
the asymptotics of the fluctuation: 



(iVex). 



C(2) 



2 r 



In 



+ 7 + l + C(2) 



(35) 



(36) 



(iii) For d = 3, the pole of the spectral Zeta function Z{j3,t) at t = 3 wins in both cases: 

(iVex)cn = C(3) (37) 

{S'N,.)en=(^]\i2). (38) 
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Of course, these results remain valid only as long as (A^'ex)cn < N. Equating {Nex)cn and 
for d = 3, e.g., one finds the large-iV condensation temperature 

for an ideal gas in a three-dimensional harmonic trap that exchanges energy, but no particles, 
with a heat bath. As expected, this result agrees with the one provided by the familiar grand 
canonical ensemble |^2|. Even more, taking into account also the next-to- leading pole, one 
obtains the improvement 

(^„)_C(3)(|^)%|C(2)(^)^ (40) 

to Eq. (PTD, implying that for Bose gases with merely a moderate number of particles the 
actual condensation temperature Tq is lowered by terms of the order N~^^^ against Tq, 



c 



C(2) 



2C(3)2/3 A^i/s 



(41) 



Even this improved canonical expression equals its grand canonical counterpart p3|-p^. 

These examples nicely illustrate the working principle of the basic integral representa- 
tions (|25| ) and (|27|): There are two opponents that place poles on the positive real axis, 
namely the spectral Zeta function Z{f3, t) on the one hand, which depends on the particular 
trap, and ({t) or ~ 1) on the other, which are entirely independent of the system. For 
both the number of excited particles and its mean-square fluctuation, the exponent of T is 
given by the location of the pole farthest to the right. Whereas the pole of ({t) and ({t — 1) 
does, naturally, not depend on the spatial dimension d, the asymptotically relevant pole of 
Z{j3,t) lies at t = d and thus moves with increasing d to the right, governing (A^'ex)cn above 
d = 1 and (5^A'"ex)cn above d = 2. 

But what about the anisotropic traps that play a major role in present experiments? 
With the ground state energy set to zero, and angular trap frequencies Ui {i = 1, . . . ,d), the 
energy levels then are 

£,yu...,ua = HuJiUi + ...+ Ud^d) = huju , z7 G INq . (42) 

The spectral Zeta function 

ziM = E (43) 

now is a Zeta function of the Barnes type (see also Ref. ||2^). Its rightmost pole is 
located at t = d, with residue 

'-'^^'"^-mi^)^' '''' 

where we have introduced the geometric mean Q of the trap frequencies. 
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(45) 



The derivation of Eq. (^^ is sketched in Appendix 0. 

The asymptotic evaluation of the canonical formulas ( pSj) and ( |27| ) now requires (3huji <^ 1 
for all i. If this condition is not met, since, for instance, one of the trapping frequencies is 
much larger than the others, one has to treat the entailing dimensional crossover effects ||2^ 
by keeping the corresponding part of Z{f3,t) as a discrete sum. In the following we will 
assume merely moderate anisotropy, so that the above inequalities are satisfied. 

For two-dimensional anisotropic harmonic traps, the computation of the canonically 
expected number of excited particles, and its fluctuation, then leads to 




C(2) 



2 r 



In 



h{uji + 0J2) ^ 



^ + ^)C(2) + /(..1,...) 



(46) 
(47) 



with 



da ae 



V 



1 -e 




(48) 



Equation (|47|) reveals a rather complicated dependence of the fluctuation {6'^Ncx)cn on 
the trap frequencies ui and 002- The comparatively simple form of the previous Eq. (|36|) for 
an isotropic trap has its technical reason in the simple expansion ( pTj ) of C(z) around its 
pole. In contrast, for two-dimensional anisotropic traps we need the analogous expansion of 
the Barnes Zeta function ( P3| ) for d = 2. The finite part of this expansion, corresponding to 
Euler's constant 7 in Eq. (^1]), now becomes a function of the frequencies ui and 002 that 
enters into the above result. Details are explained in Appendix ^, where we also show the 
identity 



/(cu,cj) = 7 + l + ln2-C(2) 



(49) 



which ensures that Eq. (^7|) reduces to the isotropic result ( p6D for ui = U2 = w. 

For any dimension d > 3, it is the pole of Z{j3, t) ait = d which determines the behaviour 
of both (iVex)cn and (5^A^ex)cn: 



(A^e.)c 



(5' AT,,), 



nn ^ 
nn 



Cid) 



Cid-i). 



(50) 
(51) 



The difference between the isotropic and the mildly anisotropic case now merely consists in 
the replacement of the frequency u by the geometric mean Q. 
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IV. A SADDLE-POINT APPROACH TO MICROCANONICAL STATISTICS 



When the ideal A^-particle Bose gas is completely isolated from its surrounding, carrying 
a total excitation energy E, one has to resort to the microcanonical framework. The micro- 
canonical counterpart of the distribution (|), that is, the probability for finding N^x out of 
the isolated particles in an excited state, is given by 



Pmc(iVex, E) 



$(iVex|^) 



Ne.<N 



(52) 



It is quite instructive to copy the previous canonical analysis as far as possible, in order to pin 
down precisely how the difference between the canonical and the microcanonical ensemble 
manifests itself. Hence, we wish to calculate the k-th moments 



N 



of this distribution 



which yield the microcanonical expectation value 



(iVcx)r 



fio{E) 



of the number of excited particles, and the corresponding mean-square fluctuation 



(53) 



(54) 



(55) 



Provided the energy E is so low that a major fraction of the particles remains in the ground 
state, i.e., if {Ne^)mc is sufficiently small as compared to A^, we can again employ the oscil- 
lator approximation (|12|), and compute the microcanonical moments fik{E) from the easily 
accessible canonical moments (^) by means of saddle-point inversions. 

To see in detail how this works, let us first carry through this program for the paradig- 
matic example of an isolated ideal Bose gas trapped by a one-dimensional harmonic oscillator 
potential. For ease of notation, we introduce the dimensionless variables a = jShuj, char- 
acterizing the inverse temperature, and n = E/{huj), corresponding to the total number 
of excitation quanta. We are thus working in the regime a <C 1, ^ 1. Again, this is 
compatible with the presence of a condensate if the particle number N is large. Writing 
fj^kin), Z{a), Si{a), and S2{a) instead of fJ.k{E), Z{j3), Si{P), and S'2(/3) (see Eqs. ( p!7| ) - 
(pjp), and defining 



Ho{a) 
Hi{a) 



1 

Si{a) 
Sfia) 



So(a) 



the inversion formula acquires the form [30 



e"""- Zia)HAa) 



1/2 



(56) 
(57) 
(58) 



(59) 



a=ak(n) 
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for = 0, 1, and 2. It is crucial to note that each moment requires its own saddle-point 
parameter ak{n), obtained by inverting the corresponding saddle-point equation 

n = --^\nZ{a) - -^\nHk{a) . (60) 
In contrast to the canonical case, we now need to evaluate the partition sum 

oo -| 

^(«) = n 1 — 7 — T (61) 

for a -C 1. This partition sum actually is a well-studied object in the theory of modular 
functions; it satisfies a fairly interesting functional equation that allows one to extract the 
desired small-a-behaviour straight away [^,^. In view of the intended transfer of the 
method to other trap types, we refrain from using this particular functional equation here, 
and resort once more to the Mellin-Barnes techniques. In this way we get 



InZ(a) = - ^ ln(l - e""^^) 

1 PT+icX) 

= — dta-*r(t)C(t)C(t+l), (62) 

ZTCt Jr—ioo 

SO that the dominant pole at t = 1 gives the approximation In Z{a) ^ ({2) /a. However, in 
order to derive a proper asymptotic formula for Z{a) we have to expand In Z{a) up to terms 
of the order 0(a°) inclusively, which necessitates to take into account also the double pole 
of the integrand (|6^) that lies at t = 0. Since the residue of this pole reads 

-C(0)lna + C'(0) = ^lna-^ln27r, (63) 
we obtain the desired approximation 

In Z(a) = ^ + -In — + 0(a) . (64) 
a 2 27r 



With Hi{a) and -^2(0) as determined by Eqs. (^) and (pif), the equations (|60D for the 
saddle-point parameters then adopt the form 



C(2) , c{k) 

n = 



a? 



(65) 



where 



c(0) = -1/2 
c(l) = -1/2 + 1 

c(2) = -1/2 + 2. (66) 



Inverting up to the required order 0(n ), one finds 

1 _ 3c(A;) 
aAn) Tc TT^ 



(67) 
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hence 



nak{n) + In Z{ak{n)) 



^ I2n If 1 ' 



(68) 



The important point to observe here is that the moment- dependent number c{k) drops out, 
so that the factors e""'=^"^Z(afc(n)) entering the inversion formula (^) become asymptotically 
equal for all k. Moreover, also the saddle-point corrections 




-1/2 



a=afe(n) 



(69) 



do not develop a significant fc-dependence, so that these parts cancel when forming the 
ratio Hence, we arrive at 



(70) 



indicating that the microcanonical expectation value for the number of excited particles be- 
comes equal to the canonical expression (^) in the asymptotic regime, where the difference 
between the saddle-point parameter ai(n) and the true inverse temperature ao(n) is negli- 
gible. Utilizing the asymptotic temperature-energy relation kBT/{huj) = v^Gn/vr obtained 
from Eq. (B^), the microcanonical counterpart to Eq. (^31) reads 



(iVex)r 



IT 



In 



TT 



+ 7 



(71) 



The calculation of the microcanonical condensate fluctuations requires more care. The 
canonical expression (^) had been an immediate consequence of the definition ([1^ and 
Eqs. (P) - (|19D, relying on the cancellation [Sf{f3) + S2{f3)] - Sf{f3) = 5*2 (/5), but now two 
different saddle-point parameters enter into the corresponding difference 



(72) 



indicating that the microcanonical fluctuation might deviate from the canonical one. How- 
ever, for the one-dimensional oscillator trap we find that Sf{a2{n)) — Sf{ai{n) is merely of 
the order ©(i/nln^n), and thus asymptotically negligible in comparison to S2{a2{n)) 
Therefore, for large n we have 



n. 



/i2(n) 
/io(n) 




= S2{a2{n)) , 



meaning 



n 



(73) 



(74) 



so that the microcanonical fluctuation of the number of excited Bose particles in a one- 
dimensional harmonic trap, considered for energies such that on the average a fraction of 
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the particles stays in the ground state, coincides asymptotically with the canonical fluctua- 
tion (g). 

This analysis can directly be translated into the language of the theory of partitions of 
integer numbers [^]. Distributing n excitation quanta among ideal Bose particles, stored in a 
one-dimensional harmonic trap, is tantamount to partitioning the number n into summands; 
a particle occupying the z/-th excited oscillator state gives a summand of magnitude u. In 



fact, utilizing Eqs. ( ]68| ) and ( |69D for computing jj,o{n) according to Eq. (|59|), one finds 



/^oH = ^-p|^.y-j , (75) 

which is just the celebrated Hardy-Ramanujan formula for the total number of unrestricted 
partitions of n |3l|, corresponding to the number of microstates that are accessible to the 



Bose particles when their common excitation energy is nhcu, provided that n does not exceed 
the particle number A^. For higher energies, the number of microstates equals the number of 
partitions that are restricted by the requirement that there be no more than summands, 
since the n quanta cannot be distributed over more than the available particles. However, 
as long as the distribution (|5^) remains sharply peaked around some value N^x ~ (A'ex)mc < 
N, even though n may be substantially larger than A^, the difference between the number 
of these restricted and that of the unrestricted partitions is insignificant, and neglecting this 
difference is nothing but the oscillator approximation ([T^) for k = 0. 

From the viewpoint of partition theory, Eq. (|7l|) gives the expected number of summands 
in a partition of n — randomly partitioning n = 1000, for instance, one expects roughly 93 
summands — and Eq. ( [74D contains the remarkable statement that the r.m.s.-fiuctuation 
of the number of parts into which n can be decomposed becomes just y/n in the asymp- 



totic limit []T6|. Higher cumulants Km{n) of the distribution which describes the number of 
summands in unrestricted partitions of n can be obtained by following the same strategy as 
outlined above for calculating K2{n) = n, leading to, e.g., 

Ksin) = — n' (76) 

and 



K,{n) = ^n\ (77) 

This indicates deviations from a Gaussian distribution, for which all cumulants higher than 
the second are zero. A general asymptotic formula for the partition moments /ifc(n) for 
arbitrary k, together with a check of this formula against exact data for k = 0, . . . ,3, can 
be found in Appendix 0. 

It is conceptually important to know that the fluctuation formula (|74D for the harmoni- 
cally trapped one- dimensional Bose gas can also be obtained without invoking the oscillator 
approximation, since the number of restricted partitions, and hence the entire distribu- 
tion (|52|) , can be well approximated with the help of an asymptotic expression due to Erdos 



and Lehner However, that approach depends on a specific asymptotic result that 

applies to the one-dimensional harmonic trap only, whereas the present method is capable 
of some generalization. 
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The preceding microcanonical analysis of the ideal Bose-Einstein condensate in a one- 
dimensional oscillator trap might appear like much ado about nothing: we have been careful 
to keep track of three slightly different saddle-point parameters, but in the end this distinc- 
tion turned out to be insignificant, and we have merely recovered the canonical results. But 
this is not true in general; the following reasoning will show that (and why) a condensate in a 
three-dimensional isotropic harmonic trapping potential behaves differently. In this case we 
again face a partition-type problem, since the total excitation energy E remains an integer 
multiple of a basic quantum hu. We can then virtually retrace the steps that have led for 
d = 1 to the microcanonical formulas ( [7I| ) and (fTil): Starting from the partition sum 



Z(a) 



u=l 



-exp(-az/)](^+i)('^+2)/2 



(78) 



and applying the Mellin-Barnes transformation, one readily finds f\ 
InZia) 



C(4) , 3C(S) , C(2)+5,„„ 



2 a » 

-iln27r + ^C'(-l) + ^C'(-2)+ 0(a) 



(79) 



which, together with the canonical expressions (37) and (|38D that now define Hi{a) and 
H2{a) via Eqs. ( |57| ) and (§^), yields the saddle-point equations 



n 



3C(4) , 3C(3) , C(2) , cik) 



+ 



+ 



+ 



(80) 



with 



c(0) 
c(l) 
c(2) 



-5/8 

-5/8 + 3 
-5/8 + 6 . 



(81) 



Again, these equations differ only to the order 0{a^^), and again the moment- dependent 
coefficient c{k) drops out when computing e"'*'=^'^''Z(afc(ri)) up to the asymptotically relevant 
terms of order O(n^): 



nak{n) + In Z{ak{n)) = 4C(4) 



n 



3/4 



.3C(4); 



+ 



+ 



C(2)-^^ 
8 C(4) 

C(3)^ C(2)C(3) 




!C(4)2 4C(4) 



^ln2vr + ^C'(-l) + ^C'(-2). 



(82) 



^This expansion had already been derived in 1951 by V.S. Nanda with the help of the Euler- 
Maclaurin summation formula, see Ref. |^5[. The Mellin-Barnes approach followed in the present 
work is much simpler, since it provides immediate access to the analytically continued Riemann 
Zeta function. 
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It is clear that this cancellation is quite general. Namely, for a d-dimensional isotropic 
harmonic trapping potential one has 



\nZ(a) 



ad+1) 



+ ... 



(83) 



and the saddle-point equations become 

dC{d+l] 



hence 



n 



n 



i/{d+i) 



-.d+l 



ak(n 



dCid+1] 



+ 



+ 



c{k) 



c{k) 



n 



-(d-i)/{d+i) 



d{d+ l)C(rf+ 1) \dCid+l] 



(84) 



(85) 



When computing nak{n) + In Z(afc(r;,)), the product nak{n) contributes a term c{k) that 
originates from the 0(a~^)-term in Eq. (P^. Otherwise, relevant /c-dependent contributions 
enter into nak{n) + lnZ(afc(n)) only via the leading terms of order 0{a~'^) that stem from 
Eq. (P^D on the one hand, and Eq. (R3|) multiplied by a on the other, summing up to 



(rf + 1) C(c/ + 1) X 



d 



n 



{d-i)/(d+iy 



X 



dC{d+l] 
-c{k) 



n 



-{d-i)/id+iy 



d{d+l)C{d 



dCid+1] 



-c{k) 



(86) 



and thus annihilating c{k). 

This little calculation, together with the inversion formula (|59|) , shows explicitly that 
the microcanonical expectation values (A'"ex)mc for rf- dimensional isotropic harmonic traps 
become asymptotically equal to their canonical counterparts: Forming the ratio (0), the 
factors e""Z(a)(— 27r9n/9a)~^/^ cancel even when evaluated at the slightly different saddle- 
point parameters ai{n) and ao(^)) so that we are left with 



(iVex), 



Si{ai{n)) . 



(87) 



The asymptotic equality of {Nc^mc and {Nex)cn then follows by observing that ai(n) becomes 
asymptotically equal to the true inverse temperature ao{n). 

The computation of the microcanonical condensate fluctuation along these lines, however, 
is a much more delicate matter. Returning to the particular example d = 3 for the sake of 
definiteness, both canonical expectation values {Ncx)cn = Si{a) and {6^Nex)cn = S2{a) are 
determined by the same simple pole of Z{f3,t) at t = 3, which means that both Si{a) and 
5*2(0) are proportional to a~^. This, in turn, implies that in contrast to the one-dimensional 
case the difference Sf{a2{n)) — Sf{ai{n)) appearing in Eq. ( [72D now is of the same order as 
S'2(a2(n)) itself, so that here the quite innocent-looking difference between the saddle-point 
equations (PPD, even though only of the order 0{a~^) and apparently hidden behind terms 
of order 0(a~^) — thus being overwhelmed much stronger than for d = 1 — must lead to an 
asymptotic difference between canonical and microcanonical condensate fluctuations: The 
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exponent of T will be the same, but the prefactors will differ. This observation forces us to 
evaluate Eq. ( |55D very carefully. We may not simply rely on the cancellation of c{k) as found 
in Eq. (p2D, but have to expand both ratios fJ,2{n) / fio{n) and /ii(n)//io(n) consistently up to 
terms of order 0{ak{n)~^) = 0{n^^^). This forces us to expand nak{n) +ln Z(ak{n)), as well 
as the saddle-point corrections, even up to terms of the order 0{n~^^'^)\ Detailed analysis 
shows that such an expansion is possible even with only the saddle-point equations (|80|) as 
input, although this is not immediately obvious. Proceeding in this manner, we find 



1 + 



33 



'3C(4)\ 

12C(4) V n ) 



3/4- 



H2{a2{n)) - 



1 + 



2 /3C(4)' 



3/4' 



C(4) V n 



Hl{a,{n)) 



S2{a,{n)) + [Sl{a,{n)) - Sl{a,{n))) + \^ (^^j 



3/4 



(88) 



The 0(n~^/'')-corrections in the square brackets arise because the inverse temperature ao(n) 
differs from the saddle-point parameters 02 (n) and ai(n); this is what causes the last term 
in the second equation. Since, moreover. 



Sl{a2{n)) - Sl{ai{n)) 




(89) 



we finally arrive at 



C(2) - T 



3C(3)^ 



C(2) - 7 



4C(4) 
3C(3)2 



4C(4) 




(90) 



This quantifies what we have anticipated: Apparently tiny differences between the three 
saddle-point parameters conspire to lower the microcanonical mean-square condensate fluc- 
tuation against the canonical result (|38|) , as a consequence of the fact that the rightmost pole 
of Z{i3,t) governs both (iVex)cn = Si{P) and ((5^A^ex)cn = S2{P). For the one-dimensional 
harmonic trap, where (A^ex)cn and (5^A^ex)cn are determined by two different poles, such an 
asymptotic difference does not exist. 

Conceptually instructive as the above calculation may be, it is also lacking elegance, 
to say the least. The reason for the appearance of cumbersome equations like Eq. (|82D or 
Eq. (|88[ ) lies in the fact that one exctracts the fluctuations from the exponentially large 
moments fik{E), taking the difference (|55D . This involves huge cancellations, as becomes 
dramatically clear already for d = 1 hj comparing the numbers listed in Tables I and V of 
Appendix 0. If one could avoid computing the microcanonical moments, and aim directly for 
the difference between canonical and microcanonical fluctuations, one should get expressions 
of a far simpler nature. The following section will show that such a strategy is actually 
feasible. 
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V. MELLIN-BARNES APPROACH TO MICROCANONICAL CONDENSATE 

FLUCTUATIONS 



We start by considering the excited-states subsystem with the fugacity z and the energy 
E as basic variables so that we have the relation N^y^ = Nex{z,E). In principle, A'^ex 
depends also on trap parameters that determine the single-particle energies, like the oscillator 
frequencies in the case of harmonic traps, but these parameters will be kept constant in the 
following. Taking the total differential. 



dz 



dz + 



then keeping the temperature T fixed, one finds 



dz 



z=l 



dz 



+ 



dE 



'dN,, 



dE 



dE , 



'dE' 
dz 



(91) 



(92) 



rJ z=i 



The left hand side equals the canonical mean-square fluctuation {S'^Nex)cn, whereas the first 
term on the r.h.s. is its microcanonical counterpart (5^A'"ex)mc- Hence, we obtain 



('5'iVex)cn - (5'iVex), 



'dN^ 
dE 



'dE' 
dz 



z=l 





z (e^lr 


2 = 1 


ksT^ (i) 


z 


2 = 1 



(93) 



Now the denominator 



kj,T^ ( ^ 1 



'dE_\ 



{6'E), 



(94) 



2 = 1 



is the canonical mean-square fluctuation of the system's energy, whereas the two partial 
derivatives in the numerator, 



'dN^ 
dT 



2 = 1 



'dE\ 
^dzl 



T 



{SN^ySE), 



(95) 



2 = 1 



both equal the canonical particle-energy correlation {6Ncx6E)cn = {NexE)cn — {Nex)cn{E)c 
Thus, we arrive at the noteworthy identity 



[{SN,x6E), 



(96) 



which expresses the difference between canonical and microcanonical condensate fluctua- 
tions in terms of quantities that can be computed entirely within the convenient canonical 
ensemble. The usefulness of this formula, first stated by Navez et al. [0, rests in the fact 
that it lends itself again to the oscillator approximation, and thus to an efficient evaluation 
by means of the Mellin-Barnes transformation: Within the oscillator approximation, the 
canonical particle-energy correlation becomes 
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d 



d 



dz 



dp 



z=l 



E 



■ / dtT{t)Z{(3,t-l)at-l) 



exp(/?e^) - 1 \exp(/?e^) - 1 

1 1 /"T+joO 

and the canonical energy fluctuation adopts the quite similar form 



(97) 



d 



^1 In /3) 



2 = 1 



E 

1 



+ 1 



\ exp(/5£:;,) - 1 \exp(/5£:^) 

1 1 /"T+joO 

-— / dtr(t)z(/3,t-2)c(t-i) 



(98) 



Hence, 



o^/;-';rdtr(t)z(/?,t-i)c(t 



^/;-'"/^dtr(t)z(/3,t-2)c(t-i: 



(99) 



Compared to the saddle-point approach in the preceding section, this formula is remark- 
ably easy to handle. Applied to the one-dimensional harmonic trap, for instance, it yields 
immediately 



1 



In 



+ 7 + 1 



Improving Eq. 



2C(2) huj 

by taking also the next-to-leading pole into account, 

2 



(100) 



C(2) 



flUJ 



1 kBT 

2 huj 



and substituting fc^T/ {huj) = ^/6n/n + 3/ (27r^) as stated by Eq. 
approximation to the microcanonical ground state fluctuation: 



n 









+ 7 + 1 





(101) 

we arrive at a refined 
(102) 



As we already know, the relative difference between canonical and microcanonical mean- 
square fluctuations vanishes asymptotically for a condensate in a one- dimensional harmonic 
trap. But still, this relative difference is of the order 0(ln^n/-yn), so that Eq. ( |102| ) is 
substantially more accurate than the previous leading-order approximation (|7^). 

This fluctuation formula ( |102D again has an interesting number-theoretical interpretation. 
As indicated in the previous section, the oscillator approximation, when applied to a one- 
dimensional harmonic trapping potential, corresponds to neglecting the difference between 
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partitions of n into no more than summands and unrestricted partitions; this remains 
exact as long as n < N. Hence, Eq. ( |102| ) provides a forteriori a fair approximation to the 
fluctuation of the number of integer summands into which the integer n can be decomposed. 
Figure 1 depicts the r.m.s-fluctuation a{n) = {d'^nojUl as approximated by Eq. (|74D (upper 
dashed line) and by Eq. ( |102|) (lower dashed line; coinciding almost with the full line), and 
compares these approximations to the exact data (full line). The latter have been computed 
numerically from the distribution (0) for the one- dimensional harmonic trap, utilizing the 
recursion relation 

min{n-A''ex,-'Vcx) 

$(iVe.|nM = E ^{k\{n-N,^)huj) (103) 

k=l 

with ^{l\lhuj) = ^{n\nhuj) = 1, assuming n < N. The agreement between the exact 
fluctuation and the improved asymptotic formula is no less than striking. It should be noted 



that, within the oscillator approximation, Eq. (102) describes the fluctuation of the ground 
state particles not only up to the "restriction temperature" Tr = {Tluj / ks) \J N / Q{2) , where 
n = N, but almost up to Tq = {hu/kB)N/ In N, where the occupation of the ground state 
becomes significant [1^ . 



The handiness of the fluctuation formula (|99D becomes fully clear when dealing with 
(i-dimensional harmonic traps, d > 2. Since isotropic harmonic traps can be considered 
as special cases, we proceed at once to anisotropic potentials, and consider d = 2 first. 
All the required technicalities have already been collected in Section ^ and Appendix 0: 
The integrand in the denominator of Eq. (|99D has its rightmost pole at t = 3 with residue 
r{3){kBT/hfiy({2); the rightmost pole of the integrand in the numerator lies at t = 4, with 
residue T{A){kBT/nnyC{3). Thus, 




implying that the difference between canonical and microcanonical fluctuations is still small, 
even if only by a logarithm, compared to the canonical fluctuations (|36D or (l47|) , respectively. 
Hence, in the asymptotic limit kBT/{hQ) 3> 1 — assuming is that large that this limit 
can be reached with a condensate — canonical and microcanonical condensate fluctuations 
still agree. But this is clearly the marginal witnessed by the fact that for d = 2 the 

poles of Z(/3,t) and C(t - 1) in Eq. fall together. 

For c? > 3, the relevant poles in the integrands of Eq. ( P5| ) are located ai t = d + 1 and 
t = d + 2, and one finds the general formula 

(|^)\ (106) 

In particular, for isotropic three-dimensional traps one recovers, but now without any sub- 
stantial effort, the previous result (|90|). More generally, for (i > 3 the condensate fluctuations 
in harmonically trapped, energetically isolated ideal Bose gases are significantly smaller than 
the corresponding fluctuations (|5lD in traps that are thermally coupled to some heat bath, 
although the exponent of T remains the same. As we have repeatedly emphasized, this 
finding is explained by the fact that for (i > 3 both {Ncx)cn and {S'^Nex)cn are determined 
by the same simple pole of the spectral Zeta function Z{i3,t). 
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VI. CONCLUSIONS 



The three central formulas explained in this work, Eqs. (P^D, (P7|), and (^), vindicate 
the assertion put forward in the Introduction: It is the location of the rightmost pole of the 
spectral Zeta function Z{(3,t) that determines the statistical properties of the condensate. 
If that pole is located at t = p, and < p < 1, then we deduce from Eq. (^) that {Nex)cn 
grows linearly with temperature (so that {no)cn decreases linearly with T) in the condensate 
regime, irrespective of the detailed properties of the trap, since in that case the pole of ((t) 
lies to the right of p. li p = 1, the poles of Z{j3,t) and ((t) fall together, so that the linear 
temperature dependence develops a logarithmic correction. If p > 1, we have {Nex)cn T^- 

The key point to be noted when discussing canonical condensate fluctuations is that the 
pole of the Riemann Zeta function ((t — 1) in Eq. ( p7D lies at t = 2, so that (5^no)cn changes 
its T-dependence at j9 = 2: If < p < 2, then {S^nQ)cn oc T^; if p = 2, there is the familiar 
logarithmic correction to this quadratic T-dependence, as expressed by Eqs. and ( ^7[ ) 
for two-dimensional harmonic traps; if p > 2, then (5^no)cra oc T^. 

The saddle-point calculations in Section |V| may be cumbersome, but they exemplify on 
an elementary level why (no)cn equals (no)m,c in the asymptotic regime, and why canonical 
and microcanonical condensate fluctuations may differ. Equation (^) summarizes the essen- 
tials for the three-dimensional oscillator trap: One needs three slightly different saddle-point 
parameters for computing the required microcanonical moments ([53|) within the oscillator 
approximation (|T2|) from their canonical counterparts {^^', these slight differences lower the 
microcanonical fluctuation against the canonical one. The elegant expression (|99|) links the 
difference {6^no)cn — (^^^o)mc again to the dominant pole of Z{P,t): If < p < 2, that 
difference has an exponent of T which is smaller than that of (5^no)cn; so that both types 
of fluctuations become asymptotically equal, but if p > 2, then the difference acquires the 
same exponent of T as (5^no)cn, so that the microcanonical condensate fluctuation remains 
lower than the canonical one even in the asymptotic regime. 

We have evaluated canonical and microcanonical condensate fluctuations explicitly for 
harmonic trapping potentials, where Z{/3, t) reduces to the familiar Riemann or Barnes-type 
Zeta functions. This may appear a bit special, but an analogous discussion is possible for 
quite arbitrary traps, if one merely exploits the connection between the residues of Z{f3,t) 
and the corresponding heat-kernel coefficients (see Ref. for a brief explanation of this 
fairly deep connection). 

The vision of letting the poles of Z{j3,t) move in the complex t-plane is not a fantasy 
restricted to the theorist's ivory tower, but may have direct experimental consequences. Con- 
tinuously deforming the trapping potential means continuously changing the trap's single- 
particle spectrum, and hence shifting p. For example, a spectrum of the type |22,n| 



e 



ui,...,ua 



SoY^d^t (106) 



with integer quantum numbers Ui and dimensionless anisotropy coefficients Cj not too differ- 
ent from unity implies p = d/s [0, so that, e.g., steepening a three-dimensional harmonic 
oscillator potential (s = 1) towards a box potential {s = 2) means lowering p from 3 to 3/2. 
During such a process, the fluctuation of a large condensate is described by an exponent 
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of T that changes as long as p remains above 2, since then ((52no)™c oc T'^/^ but remains 
constant when p decreases further; {6'^no)mc oc T^. 

There is another detail that deserves to be mentioned. From the viewpoint of parti- 
tion theory, Eq. ( f^D states that the r.m.s. -fluctuation of the number of parts into which a 
large integer n can be decomposed is approximately normal, (y{n) ~ ^/n. However, when 
characterizing condensate fluctuations, one would not do so in terms of the number of exci- 
tation quanta n, but rather in terms of the number of excited particles {Nf,^)mc Q- -But then 
Eqs. (|7I|) and yield 



{5^N,Xi'^ oc {N,,)rnc , (107) 

apart from logarithmic corrections, stating that the normal partition-theoretic fluctuation 
translates into supranormal fluctuation of the number of excited Bose particles in a one- 
dimensional harmonic trap. More generally, for traps with single-particle spectra ( |106| ) one 



obtains |T^ 



{6^N,,Y/^ oc (iVex) for < rf/s < 1 , 

(^^iVe.)^/^ oc {N.^Y'" for l<d/s<2, 

(52iV,,)V2 oc (iVe,)i/2 for 2<d/s- (108) 

both within the canonical and the microcanonical ensemble. Hence, when increasing d/s 
from the Id-harmonic oscillator value 1, the degree of supranormality is gradually lowered, 
until one arrives at normal particle number fluctuations for d/s > 2. 



Taking these insights together with those obtained in the related previous works p0|-p!6 
it seems fair to conclude that by now a classic problem in statistical mechanics, the fluctu 
ation of an ideal Bose-Einstein condensate, has been fully understood. 



APPENDIX A: RESIDUES OF BARNES-TYPE ZETA FUNCTIONS 



According to Section pl| , the spectral Zeta function for d- dimensional anisotropic har- 



monic traps adopts the Barnes form 

mt)^ i; (Ai) 

and the canonical thermodynamics of an ideal Bose gas stored in such a trap depends 
crucially on the rightmost pole of this function. In this appendix we briefly sketch the 
derivation of Eq. (@), i.e., of the residue of the rightmost pole of Z{P,t). 



The starting point is the contour integral representation of the Gamma function ||37[| 



^We choose to characterize the fluctuation in terms of {Ncx), rather than (no), since the properties 
of an ideal Bose-Einstein condensate, including its fluctuation, are independent of (no). This is 
just what is exploited in the oscillator approximation. 
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T(t) = dai-aY-^e-'' , (A2) 

^ ^ 2 sin(7rt) Jc ^ ^ ^ ^ 

where C is enclosing the positive real axis counterclockwise. With the help of this represen- 
tation we deduce 

Z(l5,t) = - /da (-a)*-^e-" V ttJ^ 

^ 2sin 7rt r t ic ^ ^ ^ ^ {[3nuuy 

2sin(7rt)r(i) ^ Jc ^ ' 



eINo/{o} 



r(i-t) 



, da (-a)*-i <^ — ^- -— -I) . (A3) 

2m Jc ^ ' lnti(l-e-"^^^0 J 

The first equality is obtained by interchanging summation and integration, then changing 
in each summand from the integration variable a to af5hLji7] the second by summing the 
resulting geometric series and utilizing the relation 

sin(7rt)r(t) = YX^- (A4) 

The poles of Z{f3,t) are featured by Eq. ( |A^ ) in a particularly transparent manner. 
Namely, the prefactor r(l — t) has simple poles at integer values t = 1, 2, 3, . . . . At these 
values the remaining contour integral may be evaluated immediately by just collecting the 
residues enclosed by C. The only possible pole contributing to the integral lies at a = 0; it 
has nonvanishing residues for t = —oo, . . . , —1, 0,1, ... ,d. Hence, the poles of Z{P, t) are 
located a.t t = 1, . . . ,d, and the residue of the rightmost pole is found to be 

d 

res Z{p, d) = i-iy-^ n iP^^iT^ res 1(1 - d) . (A5) 



i=l 



Using the identity 



we arrive directly at Eq. 



res T(-n) = ^ — ^ , n G INq , (A6) 
n! 



APPENDIX B: CANONICAL CONDENSATE FLUCTUATION FOR D = 2 

When evaluating the fluctuation formula (^) for two-dimensional harmonic traps, the 
product Z{P,t) ((t — 1) provides a double pole at t = 2. In that case the knowledge of the 
residue (0) is not enough for computing the mean-square condensate fluctuation; also the 
finite part of Z{f3,t) at t = 2 enters into the residue of the double pole. More precisely, in 
analogy to Eq. (|31|) for the Riemann Zeta function, one needs the expansion 
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for t close to 2. In this appendix we determine the function /(tui, to'2, t), and thus prove 
Eq. (13). 



Introducing a = Juji/uj2 and 6 = wu;2/^^i, we first write 



t 



^ ^ i7gINo/{o} 

valid for Be{t) > 2. Splitting the sum according to the scheme 

00 00 00 

^elNo/{o} ^^1=^ ^2=^ ^1''^^=^ 



we find the decomposition 

Z{P,t) = (^^j {C(t)(a-* + r*)+//(a;i,a;2,t)} , (B3) 

where 

00 

H{uJi,UJ2,t) = E 



1 



1 /-oo „-{a+6)a 

daa*-^- ^. (B4) 



T{t) Jo (1 - e-«°) (1 - e-''") 

This identity is obtained in a similar manner as Eq. using the familiar representation 

/•oo 

r(t) = / daa*-^e-" (B5) 



of the Gamma function. 

Now we are interested in the behaviour of Z{f3,t) as t 2, where, as we know from 
Appendix ^, it has a simple pole. How is this realized in Eq. ( |B3| ) ? Since ({t) is regular at 
t = 2, the pole is contained in the integral ( p^) . At the lower integration bound, that is, for 
a — i> 0, the integrand behaves as 1/a for t ^ 2; therefore the integral diverges at t = 2. The 
behaviour of the integral as t tends to 2 is extracted with the help of the following trick. 
For ife (t) > 2, write 

1 /"OO , s ( 1 1 1 \ 

i7(cui, t) = -—\ d« «*-ie-('^+')" - ^ + ^ 

1 /"OO , / 1 1 \ 



rm io I 1 



2 



r(t) 7o V (1 - e-"") (1 - e-^") a 

where Eq. ( P5| ) has been used. The simple pole of Z{(3,t) at t = 2 is now contained in the 
first term, since r{t — 2)/r(t) = l/[{t — l){t — 2)], and the remaining integral is finite for 
t = 2. In this way, we arrive at the expansion 
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H{ui,U2,t) 



t-2 



- 1 - In 



+ I{ui,u;2)+0{t-2) 



(B6) 



with /(ci;i,ci;2) as defined in Eq. (|^). Together with Eq. (p3|), this determines the desired 
function f{uji,uj2,t) and thereby leads to the resuh (^). 

It is quite interesting to see how the fiuctuation formula (|36D for the isotropic case is 
recovered in the limit ui = uj2 = uj. Then the integral simplifies to 



I(uj, uj) 



da ae 



-f 

Jo 



da 



-2a 



a l-e-°/ 1-e 



ae 



(B7) 



Employing now the identities 

d 



dz 



\nT(z) 



r°° (\ 1 

ln;z+ / dae""" 

7o Va 1 - e-" 



fB8) 



for the Psi function, and 



r(^) h 



da 



1-e- 



(B9) 



for the Hurwitz Zeta function, we end up with Eq. (^Qf ). This equation confirms that the 
complicated expression ( ^7]) for the canonical condensate fluctuation in a two-dimensional 
anisotropic harmonic trap indeed becomes equal to the expression (|36D in the isotropic limit. 



APPENDIX C: MOMENTS OF PARTITIONS 



The saddle-point method followed in Section |^ can be employed to derive asymptotic 
expressions for the fc-th moments Hk{n) of unrestricted partitions of integer n, for arbi- 



trary k |jT9|. Defining the symbol 

^ [Ao, Ai, . . . , Afe_i] = ^ 



Afc-i 



III 12- ■ ■ ■ tfci \ 11 / V ^' , 

where the sum extends over all partitions of k, i.e., ii + 2^2 + . . . + kik = k, we find 



(CI) 



fJ'kin) 



1 / I2n 

exp TT] 



n 




6n 



TT 



E 



In 



6n 



TT 



+ 7,C(2),2!C(3),...,(A;-1)!C(A;) 



(C2) 



For A; = 0, this expression gives the Hardy-Ramanujan formula (^); for k = 1, 2, and 3, it 
adopts the forms 
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yUi(n) ~ /io(n) 

TT 

/i2(ri) ~ /io(n) ^ 



6n 



Qn 



TT 



Qn 



TT 



TT 



In 



In 



Qn 



TT 



6r2 



+ 7 +C(2) 



+ 7 +3 In 



'Qn 



TT 



+ 7 C(2) + 2C(3) 



(C3) 

(C4) 

(C5) 



Note that our result ( |C^ ) differs for k > 2 from tlie formula stated by Richmond [jT9 
and remedies the discrepancies found by this author when comparing his formula with exact 
numerical data. In fact, the above expressions are fairly accurate; some exact values of //a:(^) 
for = to 3 are iuxtaposed in Tables I to IV to the respective asymptotic predictions. For 
completeness, exact values of the r.m.s.-fluctuation cr(n) of the number of parts occuring in 
unrestricted partitions of n are listed in Table V, together with the approximation furnished 
by Eq. (|102|) . Comparing the numbers in this table to those in Table I, one gets a vivid 
impression what it means to isolate microcanonical fluctuations from an exponentially large 
background. 
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FIGURES 



FIG. 1. R.m.s. -fluctuation cr(n) of the number of integer summands into which the integer n 
can be partitioned. The upper dashed hne is the leading approximation cr{n) ~ ^/n; the lower 
dashed line (coinciding almost with the full line) is the more accurate approximation obtained 
from the square root of Eq. ( |102| ). The full line indicates the exact values. Some numerical data 
are listed in Table V of Appendix |C|. From the viewpoint of statistical mechanics, the upper dashed 
line gives the r.m.s. -fluctuation of the number of ground state particles for an ideal Bose gas in 
a one-dimensional harmonic oscillator trap kept in contact with some heat bath, such that the 
average number n of excitation quanta does not exceed the particle number. The other two lines 
correspond to the (approximate and exact) microcanonical condensate fluctuation, that is, to the 
r.m.s. -fluctuation of the number of ground state particles when the gas is totally isolated from its 
surrounding, carrying n excitation quanta. 
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TABLES 



n 


LLnin) (exact) 


iinin] (asymptotic) 


rel. error 


50 


0.2042260 • 10^ 


0.2175905 • 10^ 


0.0654 


100 


0.1905693 • 10^ 


0.1992809 • 10^ 


0.0457 


200 


0.3972999 • lO^^ 


0.4100251 • 10^3 


0.0320 


300 


0.9253083 • lO^^ 


0.9494095 • 10^^ 


0.0260 


500 


0.2300165 • 10^2 


0.2346387 • 10^2 


0.0201 


1000 


0.2406147 • 10^2 


0.2440200 • 10^2 


0.0142 


1500 


0.1329462 • 10^° 


0.1344797 • 10^° 


0.0115 



TABLE I. Comparison of exact numbers ^o{n) of unrestricted partitions of n with the 
Hardy-Ramanujan approximation ([75|). 



n 


/ii(n) (exact) 


Hi{n) (asymptotic) 


rel. error 


50 


0.2805218 • 10^ 


0.2740428 • 10^ 


0.0231 


100 


0.4144913 • 10^° 


0.4087936 • 10^° 


0.0137 


200 


0.1357412 • 10^^ 


0.1346191 • 10^^ 


0.0083 


300 


0.4102848 • 10^^ 


0.4077577 • 10^*^ 


0.0062 


500 


0.1411488 • 1024 


0.1405470 • 1024 


0.0043 


1000 


0.2281551 • 10^^ 


0.2275624 • 10^^ 


0.0026 


1500 


0.1621438 • 10^2 


0.1618281 • 10^2 


0.0019 



TABLE IL Comparison of exact first moments Hi{n) of unrestricted partitions of n with the 
asymptotic formula (C3). 



n 


^J,2{n) (exact) 


fi2{n) (asymptotic) 


rel. error 


50 


0.4461898 • 10* 


0.4539366 • 10*^ 


0.0174 


100 


0.1027721 • 10^2 


0.1037857 • 10^2 


0.0099 


200 


0.5209742 • 10^*^ 


0.5239850 • 10^^ 


0.0058 


300 


0.2027390 • 1020 


0.2036083 • 1020 


0.0043 


500 


0.9563321 • 1025 


0.9591871 • 1025 


0.0030 


1000 


0.2361756 • 10^6 


0.2366168 • 10^6 


0.0019 


1500 


0.2146020 • 10"^^ 


0.2149100 • 10^^ 


0.0014 



TABLE III. Comparison of exact second moments ^2(^1-) of unrestricted partitions of n with 



the asymptotic formula (|C4|) . 
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n 


u^fn) (exact) 


ii-\(n) (asymptotic) 


rel. error 


50 


0.8145597 • 10^ 


0.9334154 • 10*^ 


0.1459 


100 


0.2898292 • lO^^ 


0.3173679 • 10^3 


0.0950 


200 


0.2249985 • 10^^ 


0.2390975 • 10^*^ 


0.0627 


300 


0.1120055 • 10^2 


0.1175340 • 10^2 


0.0494 


500 


0.7186145 • lO^^^ 


0.7449881 • lO^^^ 


0.0367 


1000 


0.2683336 • 10^8 


0.2749645 • 10^8 


0.0247 


1500 


0.3099702 • 10^*^ 


0.3160662 • 10^^ 


0.0197 



TABLE IV. Comparison of exact third moments ^3(n) of unrestricted partitions of n with the 
asymptotic formula (C5). 



n 


(j(n) (exact) 


cr(n) (asymptotic) 


rel. error 


50 


5.46 


5.65 


0.0349 


100 


8.14 


8.29 


0.0190 


200 


12.00 


12.12 


0.0104 


300 


15.00 


15.11 


0.0073 


500 


19.80 


19.89 


0.0047 


1000 


28.71 


28.79 


0.0026 


1500 


35.60 


35.66 


0.0018 



TABLE V. Comparison of exact r.m.s. -fluctuations (T{n) of the number of parts in unrestricted 
partitions of n with the predictions obtained by taking the square root of the asymptotic Eq. ( p.02D . 
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